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We obtain the quantum Langevin equation (QLE) of a charged quantum particle moving in a 
harmonic potential in the presence of a uniform external magnetic field and linearly coupled to a 
quantum heat bath through momentum variables. The bath is modeled as a collection of independent 
quantum harmonic oscillators. The QLE involves a random force which does not depend on the 
magnetic field, and a quantum-generalized classical Lorentz force. These features are also present 
in the QLE for the case of particle-bath coupling through coordinate variables. However, significant 
differences are also observed. For example, the mean force in the QLE is characterized by a memory 
function that depends explicitly on the magnetic field. The random force has a modified form 
with correlation and commutator different from those in the case of coordinate-coordinate coupling. 
Moreover, the coupling constants, in addition to appearing in the random force and in the mean 
force, also renormalize the inertial term and the harmonic potential term in the QLE. 
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I. INTRODUCTION 

The issue of the magnetic response of a charged quan- 
tum particle moving in a potential arises in many prob- 
lems of theoretical and experimental relevance-e.g., Lan- 
dau diamagnetism [l|, Q , quantum Hall effect pi, 0] , two- 
dimensional electronic systems [U, and others. The ad- 
ditional effect of quantum dissipation due to interaction 
with the external environment may be studied system- 
atically by employing the systernjDlus-reservoir model, 
i.e., the Caldeira-Leggett model M (also known as the 
independent-oscillator model [7H12|]). In this scheme, the 
environment is modelled as a quantum mechanical heat 
bath or reservoir comprising an infinite number of inde- 
pendent quantum harmonic oscillators with continuously 
distributed frequencies. One assumes a specific coupling 
of the dynamical variables of the oscillators to those of 
the particle. 

In the case of bilinear coupling between the particle 
coordinate and the coordinate of each bath oscillator, a 
reduced description of the particle motion is given by the 
quantum Langevin equation (QLE) satisfied by the par- 
ticle coordinate operator. In this equation, coupling to 
the bath is described by (i) an operator-valued random 
force, an d ( i i) a mean force characterized by a memory 
function [Tl|, [l3| . These forces do not depend on the mag- 
netic field whose only appearance in the QLE is through 
a quantum generalization of the classical Lorentz force. 

In this work, we consider the complementary possibil- 
ity of coupling of a quantum system to a quantum me- 
chanical heat bath through the momentum variables. Al- 
though such a scenario has been considered previously by 
many authors [6 U14Hl9l |. here we study the additional fea- 
ture of the presence of an external magnetic field. To this 
end, we consider a gauge-invariant system-plus-reservoir 



model. The system comprises a charged quantum parti- 
cle moving in a harmonic potential in the presence of a 
magnetic field. The particle is linearly coupled via the 
momentum variables to a quantum heat bath consisting 
of independent quantum harmonic oscillators. 



Here, we derive a QLE for the particle coordinate op- 
erator for the case of an external magnetic field which 
is uniform in space. The QLE is obtained by utilizing 
the well-known Hcisenberg equation of motion for evolu- 
tion of quantum operators and by effectively integrating 
out the bath degrees of freedom from the equations of 
motion. We show that similar to the case of coordinate- 
coordinate coupling, the QLE involves (i) a quantum- 
generalized Lorentz force term, and (ii) a random force 
which docs not depend on the magnetic field. This latter 
force, nevertheless, has a modified form, with symmet- 
ric correlation and unequal time commutator different 
from the corresponding results in the case of coordinate- 
coordinate coupling. Other differences include (i) the 
memory function characterizing the mean force in the 
QLE has an explicit dependence on the magnetic field, 
and (ii) the inertial term and the harmonic potential term 
in the QLE get renormalized by the coupling constants. 



The paper is organized as follows. In the next sec- 
tion, we introduce the system of study, and show that 
the system is invariant under a gauge transformation. In 
Section UlTl we derive the Heisenberg equations of motion 
for the particle and the bath oscillators. In Section ITVl 
we derive the QLE for the charged particle for the case 
of a magnetic field which is uniform in space. Finally, we 
draw our conclusions. 
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II. SYSTEM OF STUDY 



Using the Hadamard formula 



Consider a charged particle moving in a harmonic po- 
tential in the presence of an external magnetic field. The 
particle is linearly coupled through the momentum vari- 
ables to a large number N of independent quantum har- 
monic oscillators constituting a heat bath. The Hamilto- 
nian of the system is given by 
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where e, m, p, r are respectively the charge, the mass, 
the momentum operator and the coordinate operator of 
the particle, while wo is the frequency characterizing its 
motion in the harmonic potential. The jth heat-bath os- 
cillator has mass rrij, frequency Uj, coordinate operator 
q j 5 and momentum operator pj. The dimensionless pa- 
rameter gj describes the coupling between the particle 
and the jth oscillator. The speed of light in vacuum is 
denoted by c. The vector potential A = A(r) is related 
to the external magnetic field B(r) through 



B(r) = Vx A(r) 



(2) 



The relevant commutation relations for the various co- 
ordinate and momentum operators are 

[r a ,Pp] = itti>a/3, [qju,Pk(l} = ihSjkSa/3, (3) 

while all other commutators vanish. In the above equa- 
tion, Sjk denotes the Kronecker 6 function. Here, and in 
the following, Greek indices (a, f3, . . .) refer to the three 
spatial directions, while Roman indices k, . . .) repre- 
sent the heat-bath oscillators. 

Wc now show that our system of study is gauge invari- 
ant. Consider the gauge transformation 



A(r)^A'(r)=A(r) + V/(r), 



(4) 



where /(r) is an arbitrary function of coordinate r. The 
transformed Hamiltonian H' = H'(A') is given by the 
right hand side of Eq. (P} with A replaced by A' . 

Now, our system will be gauge invariant if simultane- 
ous with the transformation one can make a unitary 
transformation of the state vectors of the system, 



\m) -> W{t)) = um))-, = u- 



(5) 



such that all physical obscrvablcs remain invariant un- 
der the joint transformation. This requires that one 
should have H'(A') = UH(A)U\ where H(A) = H. In 
our case, finding such a unitary transformation is easily 
achieved with the choice 



U 



cxp (|;/(r; 
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(7) 



and the commutation relations ([3]), one can check that 
H'(A') = UH(A)U\ as required. 



III. HEISENBERG EQUATIONS OF MOTION 

In this section, we derive the Heisenberg equations of 
motion for the charged particle and the heat-bath oscil- 
lators. 



1. Charged particle 

For the charged particle, the Heisenberg equations of 
motion are 

v = r = ^-[r,H] 

in 

-A) .'/;P ' — A). 
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and 



Pa = ^r[p a ,H] 
in 



e 
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Equation 



(d a Ap)vp + vp(d a Ap) 
gives 



"0 7 o- 
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where m r is the "rcnormahzed mass" , defined as 



N 2 

L z — ' m> 

j=l 3 



Next, note that 

(v x B) Q = vpd a Ap - vpdpA a , 

and that 

(d a Ap)vfj = Vfs(d a Ap) + [d a Ap,vp] 

ih 

= vp(d a Ap) H dadpAp. 

m r 

Using Eqs. (JT^J) and (T3J) in Eq. ©, we get 

6 G tflG 
Pa = ~(V X B) Q + -VpdpA a + . 

c c 2m T c 



(11) 



(12) 



(13) 



-d a dpAp - muj r a , 
(14) 
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that is, 

p = -(v x B) + -(v.V)A + -^-V(V.A) - mw„r. 

c c 2m r c 

(15) 

Now, wc have 

A a = — [A a , H\ 
in 

ih 

= vp(dpA a ) + ^— dpdpA a , (16) 

so that 

A(r) = (v.V)A + -^V 2 A, (17) 
2m r 

which, on substituting in Eq. (|10l) . gives 

■fc N 

p = m r r + -(v.V)A + - V 2 A+> ^-^p,-. 

c 2m r c m, 

(18) 

On equating Eq. (fT5"|) with Eq. ([15]). we get 

m r f = — mus^c: H — (v x B) 

c 

+^(v(V.A)-V 2 A)-f 

2m r cV v ; / ^ m,- J 

(19) 

On noting that 

47T 

V(V.A) - V 2 A = V x (V x A) = V x B = — j, (20) 

c 

where j is the current producing the external magnetic 
field, and also the fact that in practice this current source 
lies outside the region where the charged particle moves, 
we have 

N 

m r r = -mw„r + -(v x B) - > — p 7 . (21) 

c ^— ' m, 

j=i j 

In the next subsection, we show that i>j = —mjCJjqj. 
Using this in the last equation, we get 

TV 

m r r = -mwjr H — (v x B) + gjm T uj?qj. (22) 
C j=i 

2. Heat-bath oscillators 

For the heat-bath oscillators, the equations of motion 
are 

q 777 fl - 77 



and 

p ' 777 p - 7/ 

= -mjujjqj. (24) 
Combining Eqs. {23) and (|23|). we get 

mjijj = —rrijUj^qj — gjp + ^— A, (25) 

which, on using Eqs. (fTB"j) and (fT7|) . gives 

rrijCij = —nijUijCij + g^mio^v — ^— (v x B) 

+ ^(Va-V(V.A)Y 

2?n r c V / 

(26) 

Using Eq. ([20| and the reasoning given in the sentence 
following it, we finally have 

mjcjj — m ,-J~<\ , + gjinui^r — ^— (v x B). (27) 



IV. UNIFORM B: THE QUANTUM LANGEVIN 
EQUATION 

In this section, we derive a QLE for the charged parti- 
cle interacting with the heat-bath oscillators as modelled 
by Eq. ([1]), where we now consider a magnetic field uni- 
form in space. One of the early appearances of a QLE in 
the case of coordinate-coordinate coupling between the 
particle and the heat-bath oscillators in the absence of 
magnetic field was in Ref. 0. In our case, we follow the 
program adopted in [TJ| for the derivation of the QLE. 
The essential steps are as follows. 

• Step 1: Obtain the Heisenberg equations of motion 
for the system of the charged particle coupled to 
the heat bath. Solve these equations for the bath 
variables, and substitute the solution into the equa- 
tions for the charged particle to obtain a reduced 
description of the particle motion. The solution 
will contain explicit expressions for the dynamical 
variables at time t in terms of their initial values. 

• Step 2: Make specific assumptions about the initial 
state of the system, e.g., assume that the heat bath 
was at thermal equilibrium at the initial instant 
with the bath variables distributed according to a 
canonical distribution. 

• Step 3: Show that the coordinate operator for the 
charged particle then represents a stochastic pro- 
cess in time, and satisfies a QLE. The statistical 
properties of the stochastic process arise from the 
initial canonical distribution of the heat bath. 

Step 1 has been partially carried out in Sec. IIII1 We 
now carry out the remaining part, and solve the equa- 
tions of motion for the bath variables by considering the 
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magnetic field B to be uniform in space, with compo- 



nents B x ,B y ,B z , and magnitude B = B x + By 
In this case, Eq. (|27|) has the retarded solution 

,2 n ™, ,2 



q,-(t) = q?(t) 



^ £r(f) - ^ £r(0) cos(wjt) 



d£ r{t')cos(ujj(t - 1')) 



9 j 1 



nijOjjB 



r / d£ r(i / )sin(w i (t-t / )), 
o 



where 



q£(i) = q,(0) cos(w^) + sin(wjt) 

771 jWj 

is the contribution from the initial condition, 
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mc 



(28) 



(29) 



(30) 



is the Larmor frequency of precessional motion of the 
charged particle in the magnetic field, and 





-B z 



By —B X 
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(31) 



Substituting Eq. ([28]) into Eq. ([22]). we get 



m r r + / dt'r(t')fi(t - t') + m v w%r + [i d (t)r(Q) 
Jo 

--(vxB)=F(t), (32) 



where 



AT 



(33) 



F(t)E^ 33mrW3 2 qJ(t)e(t), 
/i(t - i') = n d (t - t') + r> od (i - t'), (34) 



where /Xd, the diagonal part of function /_t, and /x ot j, the 
off-diagonal part, are given by 



N 2 2 

(t -t') = J2 m cos(u)j{t - t'))e(t - 1'), 
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(35) 



W (i - t') = V ^ ^ sin( Wj -(t - 0)e(t - O- 

(36) 



This completes step 1 of the program. 

To implement step 2, we now assume that at distant 
past, t = — oo, there was no magnetic field, the charged 
particle was held fixed at r(0), while the heat-bath oscil- 
lators were kept in weak contact with another heat bath 
at temperature T so as to be able to come to thermal 



equilibrium. Therefore, at time t = 0, the heat-bath os- 
cillators are in canonical equilibrium at temperature T 
with respect to the free oscillator Hamiltonian 



N r P 2 

y — 

^ L2m,- 

3=1 3 



(37) 



Subsequently, at a time t > 0, the particle is released and 
the magnetic field is turned on, so that further evolution 
of the system is governed by Hamiltonian ([TJ) . Note that 
this physical picture is consistent with choosing the re- 
tarded solution (f28|). The state of the system at t = 0, 
corresponding to a correlation-free preparation, is given 
by the total density matrix operator 



P t(0) = pp(0)(>W 



(38) 



where the initial density matrix operator pp(0) of the 
charged particle is given by 



pp(0)=*(r-r(0))$(p), 



(39) 



while that of the heat bath, which is in canonical equi- 
librium, is given by 



Pb = 



Z B = Tr B (e- HB/kBT ). 



(40) 



Here, k B is the Boltzmann constant. The normaliza- 
tion factor is denoted by Zb, while Tr B represents partial 
trace operation with respect to the bath variables. 

The statistical average of a heat-bath operator O with 
respect to the initial state (|38| is given by 



(O) =Tr B (Oe 



-H B /k B T 



)/Tr B (e 



-H B /k B T 



(41) 



Using known properties of quantum harmonic oscillators, 
it is straightforward to show that 

<«,•«((>)> = 0, 

(p ja (0)) =0, 



<5ia(0)g fc/3 (0)> 



(Pja(0)p fc/} (0)> 



coth ( ^-A^ )S jk 5 a f3, 



■ coth 



2k B T 
\2k B T 



5 jk 5 a 



1 



(Qja(0)Pkp(0)) = -(Pja(0)qkp(0)) = ^5 jk 5 af3 . 

(42) 

In addition, we have the Gaussian property: the statis- 
tical average of an odd number of factors of qj a (0) and 
Pj a (0) is zero, while that of an even number of factors is 
equal to the sum of products of pair averages with the 
order of the factors preserved. 

Using the results in Eq. (|42|) . one finds that the force 
operator F(t), Eq. ([33]) . has zero mean, 



(P(t)> = 0, 



(43) 
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and a symmetric correlation given by 
\(F a {t)F p (t') + F l3 {t')F a (t)) 

N _2™2. ,3 



2 e . — ' TO 

3=1 



(44) 



In addition, F(t) has the Gaussian property, which fol- 
lows from the same property of the qj(0) and Pj(0). 

Thus, the initial distribution of the heat bath oscilla- 
tors turns the force operator F(t) into an operator- valued 
random force. On using the canonical commutation rules 
©, we find that F(t) has the unequal time commutator 
given by 



N q 2 m 2 uj 3 

[F a (t),Fp(t')} = -iH5 a , J2 - J. 3 sm(Ljj(t-t')). 



3=1 

(45) 

We are now in a position to achieve Step 3 and in- 
terpret Eq. (|3"2")l with t > as a QLE for the particle 
coordinate operator r(i), which now reads 

m r r + / dt'r(t')iJ,(t-t')+m x oj 2 r + fj,d(t)r(0) 
Jo 

--(vxB)=F(t), (46) 
c 

where F(t) represents a random force with correlation 
and unequal time commutator given by Eqs. (|44|) and 
([45"T) , respectively The renormalized mass ui Y is given by 
Eq. (jlll) . The second term on the left represents a mean 
force characterized by the friction kernel or the memory 
function fj,(t). Note the appearance of the initial value 
term that depends explicitly on the initial coordinate of 
the particle and the diagonal part of the memory func- 
tion. One can absorb this term into the definition of 
the random force by defining G(t) = F(t) — ^d(i)r(0), 
and then considering the initial state (|38|) . with par- 
ticle density operator ([55]) and bath density operator 



PB 



Zb 



where the "shifted" bath Hamilto- 



nian is F|hifted 



JV 



2 m , 



\m-iUJ 2 [(\j 



9j» 



L r(0)] 



[12l |. This procedure guarantees that the redefined ran- 
dom force G(t) has the same statistical properties as 
F(t). 

We now point out some interesting features of the QLE 
(|4"6")l . which are not present in the QLE for the case of 
coordinate-coordinate coupling [l3[. These are (i) The 
coupling renormalizes the inertial mass, (ii) the harmonic 
potential term is also renormalized, (iii) the friction ker- 
nel has an off-diagonal part arising from the magnetic 
field and a diagonal part due to the harmonic potential. 
Similar to the coordinate-coordinate coupling, the mag- 
netic field appears in the QLE as a quantum-generalized 
classical Lorentz force term, and the random force in the 
QLE does not depend on the magnetic field. This lat- 
ter force, nevertheless, has a different form so that its 



symmetric correlation and unequal time commutator are 
modified from the corresponding expressions in the case 
of coordinate-coordinate coupling. 

It is interesting to see that the correlation and com- 
mutator of the random force F(t) may be related to the 
friction kernel fi{t — t'). The Laplace transform of its 
diagonal part fid(t), Eq. (|35[) . is given by 

/>oo 

Ji d (u) = dt fi(t)e iut ; Im(w) > 
Jo 

U m i Jo 



N 9 9 

9 m . 



1 



3 = 1 



UJ — UJj UJ + UJ,j 



(47) 



Using the well-known result that l/(x + i0 + ) = P{l/x)- 
iir5(x), we have 

Re /i d w + z0 + = - > -= 

/ ^ — ' m ■ 



3=1 



X (S(U) - UJj) + 8(U3 + Ldj)^j . 



(48) 



so that Eq. (|4"4l may be rewritten as 
\(F a (t)F {t>) + Fp(t>)F a {t)) 



— / cLu Rc\ii d (uj + iO )\ — 5 — 

7T J ufim 



x coth I 



(49) 



and similarly, Eq. (|3SJ) as 

[F a (t), Fp{t')} = — ^ \ duj Re[^ + i0 + )] 
™ Jo 

x sin(w(t — t')). 



CONCLUSIONS 



WgTO 



(50) 



In this work, we derived a quantum Langcvin equa- 
tion (QLE) for a charged quantum particle moving in a 
harmonic potential in the presence of a uniform external 
magnetic field and coupled linearly through the momen- 
tum variables to a collection of independent quantum 
harmonic oscillators constituting a heat bath. In this 
QLE, the magnetic field appears through a quantum- 
generalized classical Lorentz force term. The QLE in- 
volves a random force which docs not depend on the mag- 
netic field. These aspects are also present in the QLE for 
the case of particle-bath coordinate-coordinate coupling 
[l3| . However, significant differences are also observed: 
(i) The random force has a modified form with symmet- 
ric correlation and unequal time commutator different 



6 



from those in the case of coordinate-coordinate coupling, 
(ii) the inertial term and the harmonic potential term in 
the QLE get renormalized, and (iii) the memory func- 
tion characterizing the mean force in the QLE has a 
field- independent diagonal part, but also an explicit field- 
dependent off-diagonal part. 
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